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Abstract 

Possibility of spontaneous magnetization in QCD and magnetic 
properties of quark matter is discussed by evaluating the magnetic 
susceptibility within Fermi-liquid theory. The screening effects for 
gluons are taken into account to figure out the specific properties of 
the magnetic transition in gauge theories. It is shown that the static 
screening effect in terms of the Debye mass does not necessarily work 
against the magnetic instability; it promotes the instability, depending 
on the coupling constant and the number of flavors. 

Many efforts have been devoted to explore the QCD diagram on the 
density (pb)~ temperature (T) plane, which is closely related to relativistic 
heavy- ion collisions, early universe and compact stars [lj. In particular, 
it would be interesting to consider the low temperature case, where quark 
degrees of freedom are dominant over gluons and the concept of the Fermi 
sea is well denned. Various correlations are expected to be realized there 
as in condensed-matter physics; actually there have been many works about 
color superconductivity at high-density quark matter pp. 

A possibility of spontaneous spin polarization in quark matter has been 
suggested by one of the authors [2J, using the one-gluon-exchange (OGE) 
interaction. The onset mechanism is very simple but it reflects an interest- 
ing aspect of magnetic phase transition in gauge theories as QED or QCD 
[21 [31 HI [5] . The direct term gives a null contribution due to the charge (color) 
neutrality, so that the leading-order contribution comes from the exchange 
term. The Fock exchange interaction then favors the spin alignment due to 
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the Pauli principle and its infinite-ranginess gives a particular density depen- 
dence of the exchange term. If it is realized, it may give a microscopic origin 
of strong magnetic field observed in compact stars, especially magnetars [6] . 

In this work we further explore the properties of the magnetic phase 
transition in QCD by examining the magnetic properties of quark matter 
in a paramagnetic phase, by the use of Landau Fermi-liquid (FL) theory 
[3 [S] . We discuss magnetic susceptibility in detail to figure out the critical 
behavior of the magnetic phase transition in gauge theories at finite density. 
We only consider here the normal phase of quark matter as a first step for 
fully understanding the coexistence of magnetic and superconducting phases. 

It is well-known that there appear infrared (IR) singularities in the FL 
interactions in the gauge theories. To improve the IR behavior of the gauge 
interactions it is necessary to take into account the screening effects for the 
gauge field jU [7] . The inclusion of the screening effects is also required by 
the argument of the HDL resummation Since the IR singularities appear 
for quasi-particles with the co-linear momenta, the soft gluon should give 
a dominant contribution. Then the particle-hole polarization function ~ 
0(g 2 fi 2 ) should be the same order of magnitude with the energy-momentum 
of soft gluons ~ O(gfi), where g and fi are the QCD coupling constant and 
quark chemical potential, respectively. For the longitudinal gluons we can 
see the static screening by the Debye mass and the IR behavior is surely 
improved. However, there is no static screening for the transverse gluons 
and there is only the dynamic screening. Thus the IR singularities are still 
left for the interactions through the exchange of the transverse gluons. 

Applying the uniform but tiny magnetic field on quark matter, one may 
evaluate magnetic susceptibility defined as 



with the magnetization (M)/ = V~ 1 (fj Jq J d 3 xqf[(—ir x V) 2 + Ejg/}, where 
we take B//z. In the following we consider the isotropic quark matter, so 
that the first orbital contribution can be discarded and only spin degrees of 
freedom contributes to magnetization. Since gluons only carries color, we can 
consider the partial Xm f° r eacri flavor . The inverse of Xm measures curva- 
ture of the free energy F((M)f) at the origin, (xm) -1 = d 2 F/d(M) 2 \(M) f =o, 

1 We, hereafter, omit the subscript to specify the flavor for simplicity. 
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so that the divergence of xti signals the spontaneous spin polarization. The 
spin degrees of freedom of a quark with momentum k is specified by the spin 
vector a M (k) [2J, a(k) = £[z + k(z ■ k)/m{E k + m)],a°(k) = (z ■ k/m with 
= (k 2 + m 2 ) 1//2 and £ = ±1. We consider the static susceptibility for a 
constant magnetic field along z— axis. Each Fermi sea for £ = ±1 is a little 
changed in the presence of the magnetic field. Then magnetization may be 
given by the difference in the number densities 5n^_ = nk,£=i — rik,f=-i with 
the distribution function, = [exp/3(ek,£ — \i — l/2gE>(k)fi q (B) + l}~ 1 , where 
ew is the quasi-particle energy in the presence of the external magnetic field 
□ . Thus we find 

{M) = ^N C J ^g D {k)5nl (2) 

where fj, q denotes the Dirac magnetic moment, \i q = e q /2m, and the gyromag- 
netic ratio g_o(k) is defined as <?i}(k) = 2(a z — k z /E k a ) = 2[1 — k 2 z jE k [E k + 
m)]. For the weak magnetic field, dynamics among the quasi-particles around 
the Fermi surface determine the response to the external magnetic field. Ac- 
tually we find that the number-difference arises as a result of change of the 
quasi-particle energy around the Fermi surface, 

6 < = J2^^<~ e ^ 
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where ek denotes the quasi-particle energy in the absence of the magnetic 
field. Here the first term imply the explicit B dependence of the quasi- 
particle energy, while the second term is a rearrangement term peculiar to 
FL theory and represent the implicit change of the quasi-particle energy 
through the change of the distribution function, Sn^ = — n^. Thus we 
may see that some non-perturbative effects are taken into account in Eq. ([3]) 
in a self-consistent manner, even if we treat the quasi-particle interaction 
/kf,qC m a perturbative way. The quasi-particle interaction consists of two 
parts, the spin independent (/k, q ) an d independent (/k,q) terms; 

ZkCqC = /k,q + CC7k,q- ( 4 ) 



2 Generally topology of the Fermi surface is modified due to the spin- momentum cou- 
pling, but one may safely discard it in the weak-field limit. 
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From Eqs. (2), (3) magnetic susceptibility is written in terms of the FL inter- 
action, 



Xm 



\ 2 



N(T) 



1 + N(T)f° 
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with go = Jj k | =fcF dQk/4:7rgo{k)- f a is an angle-averaged Landau-Migdal 
parameter, which measures the strength of the spin dependent interaction, 



f a = 



dflu f dQ 



All 



"q fa 



47r ^k,q|| k | = | q | =fcF 



(6) 



at T = 0. 

N(T) is the effective density of states at the Fermi surface, 



N(T) 



-2N r 



d 3 k <9n k 



and is simply written as 
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in the limit of zero temperature . Eq. (jSJ) defines the Fermi velocity, which 
is given by the using the Lorentz transformation [8], 



dm 



vf 



|k|=fc F f 1 



kp_ _ N c k 2 F 

2 Jl: 
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where f( is a spin-averaged Landau-Migdal parameter defined by 
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(10) 



|k| = |q|=fc F 



with the relative angle # kq of k and q. 

Finally the magnetic susceptibility at zero temperature can be written in 
terms of the Landau-Migdal parameters, 



Xm — XPauii 



1 + 



N c k F fi 



-1 



where XPauli is the usual one for the Pauli paramagnetism, XPauii 
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= g 2 D /j, 2 q N c k F n/4:7T 2 . 



When we consider the color-symmetric forward scattering amplitude of 
the two quarks around the Fermi surface by the one gluon exchange interac- 
tion (OGE), the direct term should be vanished due to the color neutrality 
of quark matter and the Fock exchange term gives a leading contribution. 
The color-symmetric OGE interaction of quasi-particles may be written, 

1 x m m 

/kC,qC = JpZ^ / k Co,qC'6 = ^T^T^kCqO ( 12 ) 

c a,b k q 

with the invariant matrix element, 

M kc , qC = -9 2 ^ti (A a /2A a /2) M^{k, C; q, C)D, u (k - q), (13) 

c 

where M» v (k, C; q, CO = l/(4m 2 )tr [7^ + m)P(a(k))i"(6 + m)P(a(q))]; P(a) 
is the projection operator to select a state with spin vector a M [2J [3] - It has 
been well known that massless gluons often causes infrared divergences in 
the Landau parameters 0, HO] ■ 

Since the one gluon exchange interaction is a long-range force and we 
consider the small energy-momentum transfer between quasi-particles, we 
must treat the gluon propagator by taking into account the hard-dense-loop 
(HDL) resummation [9]. Thus we take into account the screening effect for 
gluons, 

D^(k -q) = PlMv) + I%D L {p) - ^ (14) 

with p = k — q, where P^ is the standard projection operator onto the 
transverse (longitudinal) mode which propagator has the form Dt(l){p) = 
(p 2 — n T (n(p)) -1 in terms of the self-energy H T ^(p). The last term repre- 
sents the gauge dependence with a parameter £. First of all we can eas- 
ily check the gauge independence of the interaction between quasi-particles. 
Since p il r p v M iiV = for the on-shell particles, any physical quantity described 
by the quasi-particle interaction is independent of the gauge parameter £ as 
should be. 

The self-energies for the transverse and longitudinal gluons are given as 
ni(po,p) = m 2 D 

t-t / \ .iru F m 2 D p k F 

iMPcP) = -« -: ■— r, u F = —, (15) 

4 |p| rLip 
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in the limit Po/|p| ~^ 0, with the Debye screening mass, m 2 D = ^2f =u d s (l/2vr 2 )5f 5 'pfk F 
[9]. Thus the longitudinal gluon is screened to have the Debye mass m^, while 
the transverse gluon receives only the dynamic screening and the Fermi-liquid 
interaction on the Fermi surface is not screened. Then the quasiparticle in- 
teraction on the Fermi surface can be written as 



2 

m 



U«c\\V\= H =k F = [-M°°D L (k - q) + M™D T (k - q)] , (16) 

with the effective coupling strength, C g = g 2 ■ 

We can see that the both Landau parameters /* , f a include the infrared 
singularities due to the absence of the static screening for the transverse 
gluons; D T (k — q) ~ — l/(k — q) 2 = —l/2k F (l — cos# k ) in this case, so 
that the logarithmic divergences appear in the Landau parameters through 
the integral over the relative angle, J df2 kq l/(l — cos# k ). 

Finally magnetic susceptibility is given as a sum of the contributions of 
the bare interaction and the static screening effect. We can see that the 
logarithmic divergences exactly cancel each other to give a finite result for 
susceptibility. 



C g N c p 



pK F L 



1 2 



(XM/XFauii) =1- °Z m(2E F + m)--(E 2 F + 4E F m-2m 2 )K\n 



2 V — p , J_- — J 
K 

(17) 

with hi = m 2 D /2k 2 F . Obviously this expression is reduced to the simple OGE 
case without screening in the limit k — > 0; one can see that the interaction 
among massless quarks gives a null contribution for the magnetic transi- 
tion. The effect of the static screening for the longitudinal gluons gives the 
contribution of g 4 ln(l / g 2 ) . In the nonrelativistic limit, it recovers the corre- 
sponding term in the RPA calculation of electron gas [U [T5] . 

To demonstrate the screening effect, we plot in Fig. 1 the magnetic suscep- 
tibility. We assume a flavor-symmetric quark matter, p u = p d = p s = ps/3, 
and take the QCD coupling constant as a c = g 2 /^n = 2.2 from the MIT bag 
model, where the value of a c is chosen so as to reproduce the mass difference 
among different spin states. Since we can see from Eq. (JTTJ) that massless 
quarks, u and d quarks, never lead to magnetic phase transition, we only 
consider the magnetic susceptibility of massive quarks, s quark with mass 
of 300 MeV in the MIT bag model. Note that all the flavors always take 
part in the Debye mass, even if we are only concerned with the magnetic 
susceptibility of s quarks. 
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Figure 1: Magnetic susceptibility at T = 0. Screening effects are shown in 
comparison with the simple OGE case: the solid curve shows the result with 
the simple OGE without screening, while the dashed and dash-dotted ones 
shows the screening effect with Nf = 1 (only s quark) and Nf = 2 + 1 (u,d,s 
quarks), respectively. 



One can see that the magnetic susceptibility diverges around kp = 1 ~ 
1.5fm _1 corresponding to order of nuclear density, where spontaneous magne- 
tization appears. The critical density for the simple OGE without screening 
is consistent with the previous one given by the energy calculation [2J. One 
may expect that the quasi-particle interaction is effectively cut off at momen- 
tum |k| = mo, which reduces the strength of the Fock exchange interaction, 
once the screening is taken into account. However, this is not necessarily 
the case in QCD. Compare the screening effects by changing the number of 
flavors (Nf); if we change Nf, the screening effect exhibits the opposite be- 
havior for Nf = 1 (only s quark) and iVV = 2 + 1 (u,d,s quarks). In the case 
of Nf = 1 the screening effect works against the magnetic phase transition 
as in QED while it favors in the case of Nf = 2 + 1. This occurs due 
to the sign change of the logarithmic term in Eq. ffTTl) at k = 2 (see Fig. 2): 
if OirNf > 2tt, k > 2 and the screening effect acts for the magnetic instabil- 
ity. □ Hence we may say that the screening effect does not necessarily work 

3 It would be interesting to consider the criterion in the large N c and Nf limits 
since a c Nf ~ 0(Nf /N c ), the screening effect changes its property, depending on the ratio 
N f /N c . 
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against the magnetic instability in QCD due to the large coupling constant 
and flavor degree of freedom. 
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Figure 2: Behavior of the logarithmic function is sketched as a function of the 
ratio, k = m 2 D /2k 2 F . Empty (full) symbols denote the case of Nf = 1(2 + 1) 
at 1.0, 1.5, 2.0fm _1 for square, circle, triangle, respectively k = 2 is the point 
characterizing the change of the qualitative behavior of susceptibility; k for 
Nf = 1 and Nf = 2 + 1 lie in the different sides. 

We have discussed the properties of the magnetic phase transition in 
QCD, by examining the magnetic susceptibility. The interesting behavior 
of the screening effect has been demonstrated, but it comes from the static 
screening for the longitudinal gluons. There is no static screening for the 
transverse gluons and the dynamic screening never contributes to the mag- 
netic phase transition at T = due to the sharpness of the Fermi surface. 
However, one may expect that the dynamic screening gives a finite contri- 
bution due to the diffuseness of the Fermi surface. Actually it is well known 
that the non- Fermi-liquid effect gives an anomalous term oc TlnT to the 
specific heat, which is responsible to the dynamic screening [12]. We will 
consider the susceptibility in the finite temperature and discuss this issue in 
another paper [14] . 
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COE "Center for the Diversity and Universality in Physics" and the Grant- 
in-Aid for Scientific Research Fund of the Ministry of Education, Culture, 
Sports, Science and Technology of Japan (13640282, 16540246). 



8 



References 

[1] T. Schafer, |hep-ph/0509068 



[2] T. Tatsumi, Phys. Lett. B489 (2000) 280. 

T. Tatsumi, E. Nakano and K. Nawa, Dark Matter, p. 39 (Nova Science 
Pub., New York, 2006). 

[3] E. Nakano, T. Maruyama and T. Tatsumi, Phys. Rev. D68 (2003) 
105001. 

T. Tatsumi, E. Nakano and T. Maruyama, Prog. Theor. Phys. Suppl. 
153 (2004) 190. 

T. Tatsumi, T. Maruyama and E. Nakano, Superdense QCD Matter and 
Compact Stars, p. 241 (Springer, 2006). 

[4] C. Herring, Exchange Interactions among Itinerant Electrons: Mag- 
netism IV (Academic press, New York, 1966) 
K. Yoshida, Theory of magnetism (Springer, Berlin, 1998). 

[5] A. Niegawa, Prog. Theor. Phys. 113 (2005) 581. 

[6] A.K. Harding and D. Lai, Rep. Prog. Phys, 69 (2006) 2631, 

[7] G. Baym and C.J. Pethick, Landau Fermi-Liquid Theory (WILEY- VCH, 
2004) 

P. Nozieres, Theory of Interacting Fermi Systems (Westview 
Press,1997). 

[8] G. Baym and S.A. Chin, Nucl. Phys. A262 (1976) 527. 

[9] J.I. Kapusta, Finite-temperature field theory (Cambridge U. Press, 
1993). 

M. Le Bellac, Thermal Field Theory (Cambridge U. Press, 1996). 
T. Schafer and F. Wilczek, Phys. Rev. D60 (1999) 114033. 

[10] T. Tatsumi, arXiv:07113349. 

[11] K. Ohnishi, M. Oka, S. Yasui, Phys. Rev. D76 (2007) 097501. 

[12] A. Gerhold, A. Ipp and A. Rebhan, Phys. Rev. D70 (2004) 105015 and 
references cited therein. 

T. Schafer and Schwenzer, Phys. Rev. D70 (2004) 054007; 114037. 



9 



[13] K.A. Brueckner and K. Sawada, Phys. Rev. 112 (1957) 328. 
B.S. Shastray, Phys. Rev. Lett. 38 (1977) 449. 

[14] K. Sato and T. Tatsumi, (2008) in preparation. 



10 



